Ion-beam driven dust ion-acoustic solitary waves in dusty plasmas 
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The nonlinear propagation of small but finite amplitude dust ion-acoustic waves (DIAWs) in an ion- 
beam driven plasma with Boltzmannian electrons, positive ions and stationary charged dust grains, is 
studied by using the standard reductive perturbation technique (RPT). It is shown that there exist two 
critical values (j c i and j C 2) of ion-beam to ion phase speed ratio (7), beyond which the beam generated 
solitons are not possible. The effects of the parameters, namely 7, the ratio of the ion-beam to plasma ion 
density (/^), the dust to ion density ratio (/x^) as well as the ion-beam to plasma ion mass ratio (/i) on 
both the amplitude and width of the stationary DIAWs are analyzed numerically, and applications of the 
results to laboratory ion-beam as well as space plasmas (e.g., auroral plasmas) are explained. 

PACS numbers: 94.05.-a; 52.35.Fp; 52.35.Sb. 



Low-temperature plasmas containing massive charged 
dust particles are frequently found in various space plasma 
environments [1] as well as in laboratory devices El 
and industrial processes [ 3 ] in the form of complex plas- 
mas. The presence of these highly massive and negatively 
charged dust particles in an electron-ion plasma is respon- 
sible for the appearance of new types of electrostatic waves 
including solitary waves, depending on whether the dust 
grains are considered to be static or mobile. These electro- 
static solitary waves have already been observed through- 
out the Earth's magnetosphere at the narrow boundaries, 
e.g., the plasma sheet boundary layer and the polar cap 
boundary layer [4]. One of these solitary waves is the dust 
ion-acoustic (DIA) wave, which is the usual ion-acoustic 
wave (I AW) modified by the presence of static dust grains. 
During the past several years, after the theoretical pre- 
diction of the existence of such DIA waves (DIAWs) by 
Shukla and Silin [5], and their experimental verification 
by Barkan et al (6), extensive works have been devoted to 
study the features of such DIAWs by many authors both 
theoretically (see e.g., [7]) and experimentally (see e.g., 

0). 

On the other hand, it has been found that sufficiently en- 
ergetic charged particles like ion-beams can significantly 
affect the propagation characteristics of solitary waves in 
plasmas (51 . In the auroral zone of the upper atmosphere, 
such types of solitary structures have been found in the 
vicinity of ion-beam regions usually having negative po- 
tentials ifTOl . The spacecraft observations in the Earth's 
plasma sheet boundary layer show the existence of both 
electrons and ions in the range of keV energy. Observations 
also indicate that both of these ions and electron beams 
can drive the broadband electrostatic waves present there 
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ifm . However, these ion beams in laboratory dusty plas- 
mas have become indispensable in the field of materials 
processing such as etching chemical vapour deposition and 
surface modification tl2j . 

A very few theoretical works on the behaviors of solitary 
waves in multi-component ion-beam plasmas have been 
done by some researchers ( see, e.g., [13]). It has been 
estimated that the presence of ion beams plays an impor- 
tant role in breaking up the solitary waves into many more 
solitons fl3l . The properties of solitary waves under the 
influence of high speed as well as slow ion beams on the 
propagation of IAWs have also been investigated both the- 
oretically and experimentally fT4l . However, the study of 
DIAWs under the influence of ion-beams, which may often 
exist in the space plasma environments (e.g., in the auro- 
ral regions), has not yet been reported in detail. Moreover, 
though the role of charged dust grains in the auroral re- 
gion has not yet been directly established, they could be 
important from theoretical view-points of space plasmas. 
Thus, the study of DIAWs under the influence of charged 
dust grains as well as ion-beams could be of interest to 
observe the ion wave oscillations in laboratory as well as 
space plasmas. This is the basic purpose of the present 
brief communication. 

We consider an unmagnetized ion-beam driven dusty 
plasma composed of positive plasma ions, positive ion- 
beams, Boltzmann distributed electrons and negatively 
charged dust grains forming only the background plasma. 
Two distributions for ions: one is the bulk, uniform cold 
ion plasma with its equilibrium flow speed equal to zero 
and other the energetic ion component, i.e., the ion beams 
having equilibrium ion-beam speed v^°\ have been con- 
sidered for the present system. The normalized set of 
basic equations describing the propagation of DIAWs is 
d t n a + d x (n a v a ) = 0, d t v a + v a d x v a = -( a d x (/) and 
= (1 + - Md) e^ + fid-jHiTib-ni, where n ai v ai 
denoting the number density and speed of a— species par- 
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tide with a — i (for plasma ions) and b (for ion beams) 
are normalized by their equilibrium values n a0 and the 
ion-sound speed c s = ^/k^T^Jm^. Here k B is the Boltz- 
mann's constant, T e is the electron temperature and m a is 
the mass of a— species particle. Also, <j) is the electrostatic 
wave potential normalized by k B T e /e, with e denoting the 
elementary charge, = (/i, 1), where // = m b /rrii 

is the ion-beam to plasma ion mass ratio. Moreover, /x d = 
Z d n d0 1 n i0 is the ratio of the equilibrium dust density (mul- 
tiplied by Zd, the number of electrons residing on the dust- 
grains) to plasma ion density and = n b0 / n i0 is the ion- 
beam to plasma ion density ratio. The space (x) and time 
(t) variables are respectively normalized by the electron 
Debye length, = yJk B T e / 47rn b0 e 2 and the inverse 
of the beam plasma frequency, uj pb = yj ' \.T\n b §e 2 / 'm b . 
At equilibrium, the overall charge neutrality condition is 
^eo + Z d n d0 = n i0 + n b0 . 

In order to derive the evolution equation for the prop- 
agation of small but finite amplitude DIAWs, we use the 
standard reductive perturbation technique (RPT) iTlSTl in 
which the independent variables are stretched as £ = 
e 1 / 2 (x — v p t) r = e 3 / 2 t . The dependent variables, on 
the other hand, can be expanded as n a = 1 + S^e^n^ , 

v a = v a0 +^T=i ejv ^ with ^o = 0, and0 = Sf =1 €^«, 
where e is a small nonzero constant measuring the weak- 
ness of the dispersion and v p is the Mach number (phase 
speed of the DIAWs normalized by the ion-sound speed, 
c s ). 

Substituting the stretched coordinates and the expres- 
sions for n a , v a and <j) into the basic equations, and 
equating the coefficients of different powers of e we get 
from the lowest order of e the expressions: — 
^/vl = (v p - v bo f , vP = ri^/vp, 

= / (v p - v b0 ) and nf ] + = 

(1 + fii — (i d ) c/)^ , together with the dispersion law v 2 = 
(fi6 2 - Hi) / (1 + \Li - n d ) 5 2 , where 5 = I-7 = 
1 — ^60 / v p is called the synchronism parameter such that 
v b o — v p for 5 = 0. Proceeding in this way we finally 
obtain the following Korteweg de Vries (KdV) equation 

d T c/) + A(j)d^(j) + Bd\(j) = 0, (1) 

where <j) = <p^\ and the nonlinear coefficient A 
and the dispersive coefficient B are given by A = 
[(1 ii d ) + 3 (/Xi - M 2 ^ 4 )] /2v p S (m - fiS 3 ) , 

B = -5 3 v 3 /2 (fn - fi5 3 ) . Note that for 5 = 0, A be- 
comes infinite and B = 0. Also, for /i^ — /x5 3 = 0, 
i.e., 5 = (^//x) 1 / 3 , both A and £> become infinite. 
In this case, DIA soliton ceases to exist, since it is 
obtained as a solution of Eq. ([]]) where A and B are 
both finite and nonzero. The stationary soliton solution 
of the KdV equation ([I]) is obtained by transforming 
the independent variables £ and r to a single new vari- 
able £ = £ — Uqt, where [/o is the constant phase 




FIG. 1: (Color online) The critical values j c i (upper panel) and 
7c2 (lower panel) of 7 are shown against \i and m . 



speed (normalized by c s ), and imposing the appropriate 
boundary conditions for localized perturbations (v/z., 
0^0, dc^/di -+ 0, d 2 ct)/di 2 -> as f -> ±00) 
as = m sec/i 2 [(£ — £/o7") /£)] . The amplitude (j) m 
(normalized by k B T e /e ) and the width D (normalized 
by \ D ) of the soliton are given by (j) m = 3U /A and 
D = y4po. 
We note that v p is real either for 7 < 1 — y/JI^JJi (when 

< 7 < 1) or for 7 > 1 — yj Hi/ n (when 7 > 1). 
Also, B < for 7 > 1 — y/fii/fi, and then the width 
of the DIA soliton becomes imaginary. So, for values of 
7 in the regime 1 — ^//^//i < 7 < 1 — ^//^//i, the ion 
beams will not be able to excite DIA solitons in our dusty 
ion-beam plasma. Moreover, B > either for 7 < 1 — 
\J or, for 7 > 1 + ^/fjbi/fjb, so that ion-beam driven 
DIA soliton excitation is possible if and only if 7 satisfies 
either 7 <mi n{l - ^W/^ 1 _ i-e. , 7 < 7 c i = 

1 - y/fJLi/fjL or, 7 > max{ l + VWMj 1 + i- e - 
7 > 7 c2 = 1 + \J \±ij \±. Thus, there exist two critical 
values of 7, namely 7 cl and 7 c2 above and below which 
the ion-beam driven DIA soliton does not exist. Both the 
critical values (one of which is less than and other is always 
greater than the unity) depend on the ion-beam to plasma 
ion density ratio as well as the ratio of their masses. 

We numerically investigate the properties of the critical 
values of 7, the phase velocity v p as well as the nonlinear 
and dispersive coefficients A and B. Figure 1 shows that 
7ci is always less than unity and it increases with decreas- 
ing the values of \±i and with increasing fi. Notice that fi > 
fii, since 7 > 0. This means that for the soliton solution to 
exist, if the ion-beam concentration increases with respect 



FIG. 2: (Color online) The phase speed v p is plotted against \ii 
and 7 for \i d — 0.8, /x = 3. 



FIG. 4: (Color online) The absolute value of the beam amplifi- 
cation rate (where 5 < 0) is shown with respect to \i and 7. The 
other parameters are as in Fig. 2. 
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FIG. 3: (Color online) Dependence of the nonlinear coefficient A 
on \ii and 7 for \i d = 0.8, /x = 3. 



to the plasma ion density, and the ion-beam mass exceeds 
that of the plasma ions, the ion-beam speed has to be very 
small compared to the phase speed. On the other hand, for 
relatively low-density of ion-beams or the smaller values 
of the ratio the beam speed may need to approach the 
phase speed for the excitation of solitons. In contrast to 
the upper panel of Fig. 1, the lower panel shows that 7 c2 
is always greater than unity, and that the ion-beam speed 
must be larger than the phase speed (when ion-beam den- 
sity increases with respect to the plasma ion concentration) 
in order to excite beam driven DIA solitons. 

Figure 2 explains the behaviors of the phase speed v p 
with respect to and 7 for = 3, ji d — 0.8. Since the 
ion-acoustic solitary wave does not couple to that of ion 
beams, one can consider a range of values of 7, which will 
enable to find v p quite easily. It is seen that v p decreases 
with decreasing the ratio 7 and for increasing the density 
ratio Hi. It also increases with increasing the mass ratio \i 
as well as with increasing the impurity parameter ji d (not 
shown in the figure), i.e. increasing the negative charge 
concentration into the dust grains. Here we have consid- 
ered an ion-beam moving in the positive direction with a 
speed greater than the critical ion-beam speed, i.e., 7 > 7 c2 
(The case of 7 < j cl is also similar) for which v p and B 
are real and finite. It is also found that at low charged dust 
impurity, v p decreases faster the larger are the ion-beam 
concentrations. The variation of the nonlinear coefficient 
A with respect to and 7 for the same parameter values 



as in Fig. 2, is shown in Fig. 3. We note that for a given 
value of v p , the amplitude of the beam-driven DIAWs de- 
pends on the various physical parameters, namely fii, 7, fi 
and ji d . From Fig. 3, it is evident that the amplitude of 
the DIA soliton is significant in the range 3 < 7 < 4. 
It decreases with increasing the ion-beam concentration. 
Physically, as the number densities of ion-beam increases 
compared to the plasma ions, the nonlinearity effect in the 
system becomes higher and higher. On the other hand, the 
mass ratio fi is also found to enhance the nonlinear coef- 
ficient A, and hence to decrease the soliton amplitude. It 
is found that that as the negative charge concentration on 
the dust grains decreases (or the nonlinear effects become 
larger) the soliton amplitude also decreases. Moreover, we 
find that the effect of the ratio 7 on the soliton amplitude is 
more pronounced at low charged dust concentration. 

Figure 4 shows that the absolute value of the beam am- 
plification rate (since 5 = 1— 7 < 0for7 > l)of 
soliton amplitude remains almost unchanged as long as the 
ion-beam concentration remains less than the plasma ion 
density, i.e., \i < 1. As the value of ji increases, the am- 
plification rate also increases, and it attains its maximum 
value at a higher ion-beam concentration. The Latter turns 
out to the increase of the amplification rate with decreas- 
ing values of 7 > 7 c2 . In the variations of the dispersive 
coefficient B with respect to the parameters as indicated 
above we find that since, B oc D, the width of the soliton, 
for a prescribed value of v p we can easily find the soliton 
widths with different plasma parameters. As for example, 
Fig. 5 shows that for a fixed ion-beam mass and constant 
charged dust concentration, the dispersive effects become 
stronger with the beam speed. As a result, the width of the 
soliton increases with increase of the ion-beam speed and 
reaches its maximum value. For relatively higher values 
of the ion-beam concentration, the rate of decrease of the 
width is comparatively high, while for lower values of the 
same, the ion-beam density does not have much effect on 
the wave dispersion, and hence the soliton width almost re- 
mains independent of it. It is found that at increasing value 
of h or for heavier ion-beams, the wave is more dispersive, 
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FIG. 5: (Color online) Dependence of the dispersive coefficient 
(B) on jjbi and 7 for the parameter values as in Fig. 2. 
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FIG. 6: Nonstationary soliton solution of the KdV equation with 
parameter values jii — 1.6, fid — 0.8, fi = 7 = 3 at r — 20 
(left panel) and r = 25 (right panel), showing the space-time 
development of the DIA solitons for an initial Gaussian positive 
pulse of the form 0=1 + exp[-(f + 2) 2 /6)]. 



and hence the width becomes higher. By reducing the neg- 
ative charge on the dust grains, one can find much higher 
soliton width than that in Fig. 5. Also, at higher ion-beam 
concentration, the width seems almost to be independent of 
Hi. It is also evident that the rate of increase of the width 
with respect to the ion-beam speed is more faster than the 
highly charged dust case. 

Numerical solution of the KdV equation ([I]) (see Fig. 6) 
shows that in the very beginning (e.g., r = 0.5), the sinu- 
soidal positive pulse propagates. As time progresses [e.g., 
r = 20, 25; see left and right panel of Fig. 6], the leading 
part of the positive pulse gets steepened due to nonlinearity 
and then as it travels more distance, the pulse breaks into a 
train of solitons due to dispersion. The small hump in front 
of the peak [left panel of Fig. 6] may be due to the reflected 
ions. The small hump appears after the peak [right panel 
of Fig. 6] and it tends to disappear after a long interval of 
time. Once the solitary peaks are generated, they propagate 
keeping their shapes unchanged due to nice balance of the 
nonlinearity and dispersion. By changing the system pa- 
rameter values one can observe different solitary peaks at 
different positions. 

To summarize, we have investigated the nonlinear prop- 
agation of small but finite amplitude DIA solitary waves in 



an ion-beam driven dusty plasma. The conditions for the 
existence of such solitary waves as well as the properties 
of the soliton amplitude and width in terms of the system 
parameters are obtained and analyzed numerically. Two 
critical values of 7 have been found beyond (above and be- 
low) which the formation of small amplitude DIA solitons 
is not possible. The predicted results could be important for 
soliton excitations in laboratory ion-beam driven plasmas 
as well as in space plasmas (e.g., in the auroral regions) 
with or without immobile charged dust grains. To con- 
clude, the DIA solitons in ion-beam driven plasmas with 
stationary charged dust grains are quite distinctive from the 
usual electron-ion plasmas, and may show experimentally 
the fascinating behaviors. Works in this direction is under- 
way, and will be communicated elsewhere. 
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